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Integral equations for thermodynamic quantities are derived in the framework of the non- 
crossing approximation (NCA). Entropy and specific heat of 4/ contribution are calculated 
without numerical differentiations of thermodynamic potential. The formulation is applied to 
systems such as PrFe4Pi2 with singlet-triplet crystalline electric field (CEF) levels. 
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The resolvent method has been developed as a povi^er- 
ful technique for the Anderson model with strong corre- 
lations.^"'' It enable us to derive dynamics as well as ther- 
modynamics. The lowest-order self-consistent approxi- 
mation is called the NCA, where the self-energy part of 
the resolvent is determined by the lowest order skeleton 
diagram. Let us first consider the simplest case where 
4/° and 4/^ configurations are relevant. Once the resol- 
vent R-yiz) for each 4/ state 7 is obtained, the partition 
function Zf of 4/ part is given by^"'' 



de e 



(1) 



where 777 (e) = —lmRy{e + iS)/n with 6 being positive 
infinitesimal is a spectral function. Physical quantities 
such as single-particle spectra and the thermodynamic 
potential — —TlnZf of 4/ contribution are derived 
with use of resolvents. 

On the other hand, entropy Sf = —dilf/dT and spe- 
cific heat Cf = TdSf/dT have been computed by nu- 
merical differentiation of rjy with respect to tempera- 
ture T. The differentiation requires one to calculate Zf 
at numerous values of temperature. To obtain reasonable 
numerical accuracy, substantial time is required in com- 
putations, especially for the second derivatives. Hence a 
method is desirable that yields the T-derivatives of Jl/ at 
a given value of temperature without numerical differen- 
tiation. In this paper, we present a method to compute 
the T-derivatives of spectral functions by solving integral 
equations, instead of taking the numerical derivatives. 

In order to establish notations for our development, 
we briefly summarize the NCA integral equations and 
their equivalents. The NCA integral equations for the 
self-energies ^^(z) of resolvents are given by'^"^ 



So(z)= I AeW{e)!{e)Y,Rp{z + e), 
Si(z) = j deW{e)[l - f{e)]Ro{z - e). 



(2) 



(3) 



where /3 denotes one of 4/^ states, /(e) is the Fermi 
distribution function, and W(e) is the product of hy- 
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bridization squared and density of states of conduction 
electron. Since the Boltzmann factor becomes singular at 
low temperatures, an auxiliary spectral function ^^(e) = 



-/3£ 



777(e) has been introduced.^ The quantity stands 
for intensity of removal of state 7. Alternatively, ^^(e) is 
interpreted as a result of the operation V, which is de- 
fined by VR{uj) = ZJ^e-'^'{-)Tr~HmR{Lu + id). Similar 
operation gives cr-yii^) — T''S-y{Lu). Then we obtain equa- 
tions equivalent to the NCA equations: 



(4) 



aoito) = / deW^(e)[l 



f{e)]^^l3ico + e), (5) 

/3 



aiiio)^ / deI^(e)/(e)^o(t^-e), 



(6) 



which gives ^^(e) without manipulation of the Boltzmann 
factor. Although linear equations (4)-(6) do not deter- 
mine the norms of ^^(w) and ay{uj), they are determined 
by the following sum-rule: 



7 



(7) 



which are obtained with use of eq. (1). 

Now we derive new integral equations for thermody- 
namics. For any energy e, Zf is represented in terms of 
spectral intensities as 



(8) 



Since 777(e) and '^7(e) are obtained by different integral 
equations, eq. (8) yields Zf. Equation (8) leads to an 
expression of entropy in terms of spectral functions of a 
resolvent 



5/= In 



T 



1 drj^ie) 1 9^7 (e) 



7,7 (e) dT dT 



(9) 



with arbitrary 7. Thus T-derivatives of 777(e) and ^7(e) 
are required. 

Performing derivatives with respect to T in eqs. (2) 
and (3) , we obtain a set of integral equations for the first 
derivatives: 



dR^iz) 
dT 



Ry{z) 



dT 



(10) 
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5Eo(^) 



^ 



dRp{z + e) 



dT J \ i-i^ I ^ Qx 



dT 



= jAeW{e)[l-f{e)]- 



dT 



+ F^{z), (11) 



(12) 



where functions Fq{z) and Fi{z) are defined as 

df{e) 



Fo{z) 



Fi{z) 



del^(e)^^i?M^ + e), (13) 



second derivatives of resolvents are obtained by follow- 
ing modifications to eqs. (10)-(12): addition of a term 
f-y{z) to the right hand side of eq. (10), and replacement 
of dRy/dT, dY,^/dT and by d^R^/dT"^, Q^S^/ST^ 
and respectively. F^{z) and /^(z) are computed with 
use of the zeroth and first derivatives of resolvent as fol- 
lows: 



■JdeW{e)^Ro{z-e), 



(14) 



^ dfie)dR0{z + e) , d^je] 
' dT 



which can be computed from resolvents. Numerical iter- 
ations of eqs. (10)-(12) yield T-derivative of resolvents. 

We now proceed to derivation of d£^^{e)/dT from inte- 
gral equations. Performing T-derivatives in eqs. (4)-(6), 

we obtain equations for d^-y{e)/dT as 



dT 

Fi{z) = -J deW{e) 

df{e) dRo{z-e) 



dT 



dT 



+ 



QX2 



d^fie 



R^iz + e) 



dT^ 



■Ro{z - e) 



dT 
duo{u)) 
dT 



da\ (w) 
dT 



= \R-y{ijj + id)\ 



.doji^ 
dT 



fdeW{e)[l-fie)]Y^ 



dT 



(15) 
(16) 



/,(z) = 2i?,(z)-^^''(^)^^-(^) 



(23) 



(24) 



(25) 



dT dT ■ 

Equations for d'^^^/dT'^ are obtained by replacement of 
d£,^/dT, da^/dT, and by d'^^^/dT'^, d'^a^/dT'^, 
and in eqs. (15)-(17), respectively. Auxiliary func- 
tions G-~f{u)) and g-y{uj) are defined by 



= j deW{e)f{ef-^^^ + G,{u;), (17) Go(a;) = -J deW{e) 



where auxiliary functions are defined by 
'dRy{LJ + i6) 



Goioo) 



2Re 



deWie) 



dT 

df{e) 



Rj{Lo + i5y 



dT 



IIC/3(^ + e), 

/3 



Gi(c.)= / deW{e)^^o{uJ-e). 



(18) 
(19) 

(20) 



^df{e)d^0{LO + e) , d^fje) 
' dT 



Gi(w) = j deWie) 

df{e)d^o{u^-e) , d^fie] 



Since eqs. (15)-(17) coincide with eqs. (4)-(6) in the case 
of vanishing functions G-y{w) and g-y{w), general solu- 
tions are given in terms of an arbitrary parameter c as 
d£,^{(jj)/dT + ci^{io) and d(j^{oj)/dT + c(j^{uj). The phys- 
ical value for c is fixed by the following condition: 

d^^{u) 



dT dT 

d\R^{io + i5)\^ da^{uj) 



■6(w - e) 



+ 



dT 

d-'\R^{u + i5)\'' 

QJ<2 



dT 
a-y{u>). 



(26) 



(27) 



(28) 



Jd.^. 



dT 



0, 



(21) 



which follows from the sum rule given by eq. (7). 

Equations for the second derivatives are derived in a 
manner similar to the first derivatives. Specific heat C/ of 
4/ contribution is represented in terms of spectral func- 
tions as 



Gf = 2T 



1 dr,^{e) 1 d^^ie) 



ri^ie) dT C^(e) dT 



dT]y{e) 



+ T^ 



r7^(e)2 V dT 

1 d'rj^ie) 



^7(e)' 



1 d%{e) 



dT 



(22) 



_V^{e) dT^ ^^(e) dT^ 

We note that C f here is the specific heat with the chem- 
ical potential fixed. A set of integral equations for the 



The particular solution with a proper coefficient corre- 
sponding to c used for the first derivatives can be ex- 
tracted in a similar way. 

We apply the above formation to the CEF singlet- 
triplet system taking PrFe4Pi2 as a target. Considering 
hybridization between 4/^ and conduction electrons of 
the molecular orbital with a„ symmetry, which only has 
spin degrees of freedom, exchange interaction is written 
as^ 

H^t = AcEF-Si • ^2 + (JiSi + J2S2) ■ Sc, (29) 

where the CEF singlet-triplet states are represented by 
dimer of two pseudo-spins 5'i and S2J and Acef is 
the CEF splitting. Wc have derived integral equations 
for resolvent and dynamical quantities in the NCA.^ Its 
application to thermodynamics in the above framework 
is more complicated but straightforward. We shall show 
only the numerical results for entropy and specific heat. 
It is known that the NCA is justified in the large n 
limit, where n is scattering channel. In this model, how- 
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Fig. 1. Entropy Sf of 4/ contribution for several values of the 
CEF splitting in the NCA and for no CEF splitting case in the 
NRG. The Kondo temperature defined by eq. (30), t4^<='»1'°s) ^ 
and that estimated from specific heat, T^'^^\ are shown. 



ever, conduction electron only has spin degeneracy, i.e., 

n = 2, whereas local 4/^ configuration has four levels. 
Hence derivation of the accurate energy scale including 
the Kondo effect cannot be expected in the NCA. If one 
can derive the acctiratc scale by another method such as 
the numerical renormalization group (NRG), however, 
one may hope that the temperature dependence can be 
reasonably understood by rcscaling the NCA results. 

Figure 1 shows temperature dependence of entropy S / 
for several values of the CEF splitting Aqef- We have 
taken the same values of coupling constants Jipc = 0.2 
and J2PC = and band width D = lO^K as those in ref. 
8. Without the CEF splitting, entropy tends to roughly 
In 2 at low temperatures. This is because only the en- 
tropy of the pseudo-spin Si disappears due to the Kondo 
screening , and the pseudo-spin ^2 remains free. On the 
other hand, there is no residual entropy in the cases of 
AcEF = 20, 40K. This means that the fixed point of the 
CEF singlet is correctly reproduced in these cases. 

Figure 2 shows specific heat C/ of 4/ contribution in 
the same parameters as those in Fig. 1. The peak in the 
case of AcEF = OK is due to the Kondo effect and that 
with AcEF = 20, 40K is due to the CEF excitations. We 
notice that the CEF splittings are about half of bare ones 
because of renormalization effect of the splitting. It is 
consistent with the numerical results for the temperature 
dependence of resistivity and the dynamical magnetic 
susceptibility.^ 

The result of specific heat C/ can be utilized for an es- 
timation of the Kondo temperature Tk- It is known that 
Cf become maximum at about Tk/3. From Cf without 



the CEF splitting (Fig. 2), we estimate as 



(NCA) 



9K, 



which is larger than the value estimated from the mag- 
netic relaxation rate.^ On the other hand, the second 
order scaling theory gives an analytical expression of the 
Kondo temperature as 



Tk = -DVM;exp(-l/Jpc), 



(30) 



which gives j'^'^'^'^'^s) ^ 30K. That energy scale is consis- 



Fig. 2. Specific heat Cf of 4/ contribution and Cf/T (inset) for 
several values of the CEF splitting computed in the NCA. 



tent with the NRG result for entropy (Fig. 1). Although 
the NCA fails to give the accurate energy scale because of 
the small number n = 2 of scattering channels, it seems 
that the physical quantities are reasonably derived ex- 
cept for very low temperatures provided they are scaled 
by the alternative energy scale . 

We now discuss application of the present theory to 
Pr skutterudites. For example, PrFe4Pi2 shows a non- 
magnetic order at low temperatures.^ Heavy fermion 
states are observed in a case where the order is sup- 
pressed by applying magnetic field^ or another case di- 
luting Pr concentration.^" Our results obtained with 
the impurity model are applicable to the latter case. 
Applying the quartet model proposed in ref. 11 to 
Fe4Pi2, we attribute the large value of Cf/T to 
the Kondo effect as shown in the inset of Fig. 2. Namely 
the degeneracy of the CEF singlet-triplet levels and the 
Kondo screening within the quartet states are relevant 
to the heavy mass in Pr2;Lai_2;Fe4Pi2. 

In conclusion, wc have presented the new set of inte- 
gral equations which give the specific heat and entropy 
in the NCA without numerical differentiation. It turns 
out the Kondo energy in the case of the singlet-triplet 
CEF system is underestimated in the NCA. With proper 
rescaling of the energy, however, the temperature depen- 
dence of thermodynamic quantities behaves reasonably 
well in the wide temperature range. 



1) Y. Kuramoto: Z. Phys. B 53 (1983) 37. 

2) N. Grewe: Z. Phys. B 53 (1983) 271. 

3) H. Keiter and G. Morandi: Phys. Rep. 109 (1984) 227. 

4) N. E. Bickers: Rev. Mod. Phys. 59 (1987) 845. 

5) Y. Kuramoto and H. Kojima: Z. Phys. B 57 (1984) 95. 

6) J. Otsuki, H. Kusunose and Y. Kuramoto: J. Phys. Soc. Jpn. 
74 (2005) 200. 

7) R. Shiina: J. Phys. Soc. Jpn. 73 (2004) 2257; R. Shiina and 
Y. Aoki: J. Phys. Soc. Jpn. 73 (2004) 541. 

8) J. Otsuki, H. Kusunose and Y. Kuramoto: J. Phys. Soc. Jpn. 
74 (2005) 2082. 

9) Y. Aoki, T. Namiki, T. D. Matsuda, K. Abe, H. Sugawara and 
H. Sato: Phys. Rev. B 65 (2002) 064446. 

10) Y. Aoki, T. Namiki, T. D. Matsuda, H. Sugawara and H. Sato: 
J. Phys. Chem. Solid 63 (2002) 1201. 

11) A. Kiss and Y. Kuramoto: J. Phys. Soc. Jpn. 74 (2005) 2530. 



